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ABSTRACT 

Context. B-type stars are known to rotate at various velocities, including very fast rotators near the critical velocity as the Be stars. 
Aims. In this paper, we provide stellar models covering the mass range between 1.7 to 15 Mq, which includes the typical mass of 
known Be stars, at Z = 0.014, 0.006, and 0.002 and for an extended range of initial velocities on the zero-age main sequence. 
Methods. We used the Geneva stellar-evolution code, including the effects of sheUular rotation, with a numerical treatment that has 
been improved so the code can precisely track the variation in the angular momentum content of the star as it changes under the 
influence of radiative winds and/or mechanical mass loss. 

Results. We discuss the impact of the initial rotation rate on the tracks in the Hertzsprung-Russell diagram, the main-sequence (MS) 
lifetimes, the evolution of the surface rotation and abundances, as well as on the ejected masses of various isotopes. Among the new 
results obtained from the present grid we find that 1) fast-rotating stars with initial masses around 1.7 Mq present at the beginning 
of the core hydrogen-burning phase quite small convective cores with respect to their slowly rotating counterparts. This fact may be 
interesting to keep in mind in the framework of the asteroseismic studies of such stars. 2) The contrast between the core and surface 
angular velocity is higher in slower rotating stars. The values presently obtained are in agreement with the very few values obtained 
for B-type stars from asteroseismology. 3) At Z = 0.002, the stars in the mass range of 1.7 to 3 M© with a mean velocity on the MS 
of the order of 150 km • s~' show N/H enhancement superior to 0.2 dex at mid-MS, and superior to 0.4 dex at the end of the MS phase. 
At solar metallicity the corresponding values are below 0.2 dex at any time in the MS. 

Conclusions. An extended database of stellar models containing 270 evolutionary tracks is provided to the community. 
Key words, stars: general - stars: evolution - stars: rotation - stars: Be - stars: mass-loss 



1. Introduction 

In the framework of the new grids of stellar models that the 
Geneva stellar evolution group has recently made available to 
the community (Ekstrom et al. 2012, hereafter Grids I), we pro- 
pose a more detailed study of some subsamples of the full grids, 
offering a reduced mass range, but a more extended coverage in 
initial rotation rates. 

An interesting subsample is the mass domain between 1.7 to 
15 Mq, which corresponds roughly to the domain of the early 
A- and B-type stars. This range of mass covers the transition 
from the massive and energetic O-stars, dominated by their stel- 
lar winds, to the less massive A stars with negligible winds. 
Then, in the B-type star range, the interplay between stellar 
winds and rotation will together affect the evolution of the star 
and the injection of mass and energy to the circumstellar envi- 
ronment, and also determine the post-main-sequence fate of the 
star. Moreover, this range contains the interesting case of Be- 
type stars. Since rapid rotation and reaching the critical velocity 
is thought at least partially to be involved to explain the Be-star 
phenomenon (Martayan et al. 2007, 2010), we provide in this 
grid an extended coverage of initial rotations up to the fastest 
cases. 



Previous studies have already explored the evolution of B 
stars towards the critical velocity as a function of the metallic- 
ity (Ekstrom et al. 2008). However, the Geneva stellar-evolution 
code at that time could not follow the model through a critical- 
rotation phase while keeping precise track of the stellar angular- 
momentum content, and computation of the models was stopped 
when they reached the beginning of that stage. Several recent im- 
provements in the Geneva code now make it possible to compute 
models of critically rotating stars, and to give a theoretical esti- 
mation of the equatorial mass loss that such stars undergo during 
this phase. In that context, we performed a new analysis of stars 
in the range of 1.7 to 15 Mq at various metallicities: Z - 0.014 
(solar), Z = 0.006 (LMC), and Z = 0.002 (SMC), with initial 
rotation rates on the zero-age main sequence (ZAMS) between 
< Q/Qcrit < 0.95. 

The present theoretical database will allow synthetic clus- 
ters to be built with the inclusion of any distribution of initial 
rotational velocities. In a forthcoming paper, we shall discuss 
our population-synthesis code and compare the results with ob- 
served stellar clusters and field star populations. 

This paper is organised as follows. In Sect. 2, we briefly re- 
call the effect of the rotation on the stellar surface properties. 

We explain the numerical method used to estimate the equato- 
rial mass loss for critically-rotating stars. In Sect. 3, we briefly 
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describe the parameters and the physical ingredients used to per- 
form our numerical simulations. The results are presented in 
Sect. 4. Some aspects of the advanced phases are discussed in 
Sect. 5. Finally, Sect. 6 contains our conclusions. 

A complete discussion concerning our faster rotating models 
and comparison with observations of Be-type stars is the subject 
of the second paper of this series (Granada et al. 2012, hereafter 
Paper II). 

2. Angular momentum losses due to stellar winds 
and equatorial mechanical mass losses: 
numerical method 

2.1. Recall of the effects of rotation on the properties of the 
stellar surface 

Rotation strongly affects the stellar surface in several ways. First 
of all, the shape of the surface itself is modified: the equatorial 
radius becoming larger than the polar one under the action of 
the centrifugal force. The shape of the star depends on several 
parameters (ratio to the critical velocity, latitudinal differential 
rotation, etc., see Zorec et al. 2011). Here, given the shellular 
hypothesis used in our stellar-evolution code (see Grids I, for 
details) the shape of the star is accounted for in the frame of the 
Roche model, as it allows for a quick numerical treatment. In this 
context, the only pertinent parameter for describing the stellar 
surface is the ratio of the actual angular velocity to the critical 
one 0} = Q/i^crit (see e.g. Georgy et al. 2011). The parameter 
Qcrit is defined as 

/ GM 

\ e, cnt 

with /?e, crit the equatorial radius when the star rotates at the an- 
gular velocity Qcrit- In the framework of this model, the maxi- 
mum ratio of the equatorial-to-polar radius (when O = Qcrit) is 
1.5. This is observationaUy supported by the observations of the 
fast-rotating star Achemar (Vinicius et al. 2006; Carciofi et al. 
2008). 

According to Maeder & Meynet (2000), the total accelera- 
tion at the surface of a rotating star, including gravity, radiation- 
pressure acceleration and centrifugal acceleration can be ex- 
pressed as (vectors are in boldface) 

gtot=geff(l-r(O,0)), (2) 

with geff = ggrav+gcen the sum of the gravitational and centrifugal 
accelerations, and r(£l,0) the local Eddington factor, account- 
ing for the effects of the rotation. Following Maeder & Meynet 
(2000), we define a velocity that vanishes the expression (2) as 
a critical velocity. The first critical velocity is given by geff = 0, 
and can be expressed as 

/2 GM 

with /?p_ crit being the polar radius when the star reaches Qcrft- 
The second critical velocity is reached when 1 - r(Q, 6) = 0. It is 
only possible for stars with a very high Eddington factor (FEdd ^ 
0.639). In this work, all the considered stellar models remain 
far from this value, and the only pertinent critical velocity to 
consider is the first one. In the following text, we systematically 
use the term "critical velocity" for "first critical velocity". 

Rotation not only affects the shape of the surface, but also 
produces latitudinal variations in the radiative flux, effective 



temperature, and mass flux, producing an anisotropic stellar 
wind (see e.g. Maeder & Meynet 2000; Maeder 2002). The ef- 
fects of such a wind on the evolution of stars remaining far from 
the critical velocity (which is the case for most of stars) are 
moderate (Georgy et al. 201 1), and can usually be neglected. In 
this paper, however, some of our models evolve very close to or 
even at the critical velocity. The effects of anisotropic winds are 
thus stronger (more than 10% in terms of angular-momentum 
loss compared to models where such effects are neglected, see 
Georgy et al. 2011), and are accounted for in this work. Also, 
wind anisotropics may produce, at least for some periods, non- 
spherical shape bubbles around rapidly-rotating stars, a feature 
that may be interesting for comparisons with observations. 

2.2. Computation of the equatorial mechanical mass loss 

During the stellar lifetime, both the surface angular velocity and 
the critical one evolve, making the rotation rate to = nsurf/Qait 
change as a function of time. It is thus possible for a star to reach 
the critical velocity during its life, even if it started on the ZAMS 
with (jjini < 1 . Once this limit is attained, the effective gravity at 
the equator of the star is nil. We thus expect a strong enhance- 
ment of the mass loss in the equatorial region, hereafter called 
"mechanical mass loss", which removes the overcritical layers 
and maintains the surface at the critical velocity or slightly be- 
low. The exact details of this process is to date not well known, 
and would need complex simulations that combine hydrodynam- 
ics and radiative transfer. It is currently not possible to perform 
such computations coupled with a stellar-evolution code for the 
full duration of the stellar life. 

To account for this mechanical mass loss in our stellar- 
evolution code, we proceed as follows^ : we estimate the angular- 
momentum amount AXrad that the model loses through radiative 
stellar winds during a time step, on the basis of the estimated 
surface quantities (radius, mass-loss rate, etc.). If the effects of 
the anisotropic stellar winds are not accounted for, we have 

2 , 

AXrad = 2^^radi^surf<urf' (4) 

with AMrad the total amount of mass lost during the time step due 
to radiative stellar winds, and rsurt the radius of the star. (This 
mean radius is defined as the radius that would have a spherical 
star having the same luminosity and effective temperature.) If the 
wind anisotropy is accounted for, this expression is more com- 
plex and needs a numerical integration over the stellar surface 
(see Georgy et al. 201 1). 

Knowing the initial characteristics of the model at the begin- 
ning of the time step, we have to estimate whether the radiative 
mass loss is sufficient to keep the star below the critical velocity 
at the end of the time step. Moreover, because of numerical diffi- 
culties, it is not possible to maintain the model exactly at the crit- 
ical velocity. We thus define a maximal ratio comax = ""n*"""'" • 
In our computations, it is set to Wmax - 0.99. 

At any time, the structure of our numerical model is com- 
posed of an envelope, where convection is non-adiabatic and ion- 
isation not complete, and of an interior zone where convection is 
adiabatic and ionisation complete. The envelope, which encom- 
passes a very small fraction of the total mass, typically one ten- 
thousandth, is assumed to rotate with the same angular velocity 



' This approach is similar to the one described in Grids I ensuring the 
conservation of angular momentum. 
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as the uppermost layers of the interior zone^. The mass removed 
by the mass loss will be a fraction of the mass of the envelope. 
The removal of these layers will remove angular momentum and 
a new angular velocity distribution inside the star will be built up 
so that the total angular momentum of the star will be decreased 
by the exact amount lost by the star. Since angular momentum 
can be transported over a given distance in a time step, the zone 
where Q is modified covers a zone with a depth approximately 
equal to the radial component of the meridional current multi- 
pUed by the time step. We thus have to compute the correcting 
factors of Q in that zone so that the process keeps precise track of 
the angular momentum. In the following, we denote A^con- as the 
number of layers that are affected by the mass loss. To estimate 
the correcting factors, called q below, we suppose that: 

- the changes in the structure of the star are small enough dur- 
ing the time step considered that they can be neglected for 
computing the q values; 

- the angular momentum removed by the radiative stellar wind 
is equally removed from the envelope and the first A^con lay- 
ers (see Georgy 2010 and Grids 1 for more details). That 
means that q is constant over the whole zone where Q will 
be affected by the mass loss, or that Qi, new = ^^i, ini (1 + q), 
with / the number of the shell considered; 

- no other angular momentum transport mechanisms act dur- 
ing the time step. 

The initial angular momentum content of the A^con- layers is 
-CtoLini = -Ce.ini + H^T-C/.ini, wherc the subscript "e" is for 
the envelope, and / concerns the /-th layers below the surface. 
Assuming that the mass carried away by the stellar winds is re- 
moved in the envelope, we can thus write the estimated total 
amount of angular momentum in the Ncoa layers after the mass 
loss occurred as 

(A \ ycorr 

1 ^^\ + {\+q)y£i^,^, (5) 
We, ini/ .^j 

with 

^e, ini the mass of the envelope at the beginning of the time 
step. Finally, since we must have Xtot, fln = Xtot, ini - AXrad> we 
obtain 



<1 = 



-e, im 



.-A 



(6) 



Knowing q, it is now possible to estimate the surface angu- 
lar velocity Qsurf, new at the end of the time step. If Qsurf, new < 
Qmax allowed^ the model does not need to lose more angular mo- 
mentum and the angular momentum losses estimate is complete. 

If i^surf, new > ^^max aUowed> the radiative mass loss is not strong 
enough to maintain the star below the critical velocity, and other 
expressions for Xtot, fin and q have to be obtained. In contrast 
to previous case, where the mass loss rate imposes an angular 
momentum loss, we need here to remove an excess of angular 
momentum, and this will determine an additional mass loss. 

We assume that this additional mass, AMmecii> is lost mechan- 
ically at the equator, and rotates at the same angular velocity as 
the surface. It thus removes the following angular momentum: 



^-£mecli ~ AA/niech^surf'" eq' 



(7) 



^ Note that during the red (super)giant phase, the envelope is much 
bigger than during the MS. When the star crosses the Hcrtzsprung- 
Russcll diagram and becomes red, the mass-coordinate of the base of 
the envelope is progressively decreased, up to M^IM^i = 0.98. 



with Teq the equatorial radius. The mechanical mass loss has 
to bring the angular velocity of the stellar surface down to 

fJmax allowed- We thuS define the factor qum as Qsurf, new = 

^max allowed - f^surf,ini(l +q\im), and the new angular velocity 
of each layer where the correction is appHed becomes 



■ i^i,ini(l +?lim) = ^i. 



a 



max allowed 



a 



(8) 



'surf, : 



As before, we assume that the efi'ect of the mechanical mass 
loss will afl'ect A^con- layers below the surface. The angular mo- 
mentum of the star that is needed to have a sub-critically rotating 
surface at the end of the time step is estimated to be 



■Ctot, Mm — -Ce, ini (1 + ^lim) I 1 ~ 



^•COTT 

+ (1 +9lim)^X;, 



AM, 



rad + AMniech \ 
^c, ini / 



(9) 



1=1 



Knowing that Xtot, lim = Xtot, ini - AXrad - AXmech, and using 
relations (7) and (8), the mechanical mass lost during the time 
step AMniech, which is the only unknown of the equations, can 
now be defined as 



AM, 



' nie.ini fisuif,! 



-AX, 



rad 



mech ■ 



^^surf, ini'^ eq 



ini ^surf , ii 



(10) 

providing an estimation on how much mass should be removed 
from the equator of the star in order to keep the star below the 
critical velocity. This mass will probably form an equatorial cir- 
cumstellar disc. In this paper, we assume that this mass is lost, 
and is not re-accreted on the star. 

As in Grids I, A'corr is set to 200 in this work. This value en- 
sures that during a characteristic time step, the difiiision and the 
advection of the angular momentum afifect at least these layers. 
To avoid a strong discontinuity at the edge of the zone where the 
correction is appUed, we limited the time step in order to Umit 
the value of q to 0.005 (0.01 for models rotating with oj < 0.2). 
Moreover, during the red (super)giant phase, during which a 
large surface convective zone develops, the correction is applied 
only to the whole convective zone, as soon as it stretches over 
more than 50 layers. 

The above method for computing the variation in Q due 
to the mass loss, while allowing the conservation of the angu- 
lar momentum, is probably quite schematic, especially for what 
concerns the equatorial mass loss when the surface velocity be- 
comes overcritical. Among the limitations, we can mention that 
the process that ultimately pushes the matter out is absent. This 
process can be momentum input by radiation and/or pulsations 
and/or convection. To fully describe such processes, other tools 
than the present stellar-evolution code are necessary. The dis- 
cussion about to what extent the results of the present models 
depend on the computational time step and on the value of Wmax 
is presented in Paper II. 



3. Physics of the models 

3.1. General context 

The physical ingredients included in our models are the same as 
in Grids I, except for the following points: 
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Fig. 1. HRD for tiie 7 Mq models with uj^ = flM/Sicnt = (black), 0.1 (blue), 0.3 (cyan), 0.5 (green), 0.6 (yellow), 0.7 (orange), 0.8 (red), 0.9 
(magenta), and 0.95 (purple). Models at Z = 0.014 {left), Z = 0.006 {centre), and Z = 0.002 {right). 



- Contrary to Grids I, some of the models presented in the 
present paper evolve for a significant time near the critical 
velocity. To follow the evolution of the surface velocities as 
best as possible, as well as the total angular momentum con- 
tent, the effect of anisotropic winds are accounted for as in 
Georgy et al. (2011). 

- Some of our models reach the critical velocity during their 
evolution. In that case, the mechanical equatorial mass loss 
is accounted for as described above. 

Our grid of models contains ten different masses (1.7, 2.0, 
2.5, 3.0, 4.0, 5.0, 7.0, 9.0, 12.0 and 15.0 M©), with nine differ- 
ent initial rotation rates (w = 0, 0.1, 0.3, 0.5, 0.6, 0.7, 0.8, 0.9, 
0.95), and this for three different metallicities: Z = 0.014 (so- 
lar metallicity), Z = 0.006 (Large Magellanic Cloud metallicity) 
and Z = 0.002 (Small Magellanic Cloud metallicity). This grid 
thus covers a mass range centred on the B-type stars^. The mod- 
els are evolved up to the helium flash (Mini < 2 Mq), the early 
asymptotic giant branch (2.5 Mq < Mini S 9Mq) or the end of 
central carbon burning (Mini ^ 9 Mq). 

For the non-solar metallicities, the initial abundances are set 
as follows 

- The hydrogen and helium abundances are computed assum- 
ing that they vary linearly with respect to the total metal con- 
tent Z from the values given by the Big Bang nucleosynthesis 
to the present values. We thus have: 



Y(Z) 



AY 
Yp + — Z, 
AZ 



(11) 



with Yp - 0.2484 the primordial helium abundance (Cyburt 
et al. 2003), and AF/AZ = 1.257. The hydrogen mass frac- 
tion is given hy X = I - Y - Z. 
- The relative abundances of all other chemical species are the 
solar ones (the abundances are only scaled on a lower total 
metal content). 

3.2. Note on the diffusion coefficients and tiie caiibrations 

Stellar evolution computations with one-dimensional codes such 
as those presented in this paper account for physical processes 



^ As in Paper II, we consider that a star is a B-type star if it has an 
effective temperature such as 10000 K < T^n < 30000 K 



that are not unidimensional (convection, rotation, turbulent mix- 
ing, etc.). One accounts for these effects through simplified the- 
ories, which have various free parameters that have to be cali- 
brated on observations: overshoot parameter, mixing length, and 
shear mixing efficiency. Moreover, the implementation of the ro- 
tation includes two diffusion coefficients: the horizontal one 
and the shear one Z)shear- In the literature, there are three differ- 
ent expressions for (Zahn 1992; Maeder 2003; Mathis et al. 
2004) and two for D^heai (Maeder 1997; Talon & Zahn 1997). 
Combining them allows for six different implementations of the 
rotation, and there is to date no reason to prefer any one of them 
on the basis of theoretical considerations. 

The calibration of the overshoot parameter and of the mixing 
length is detailed in Grids I. In the framework of this paper, we 
add some complementary details on the choice and the calibra- 
tion of the diffusion coefficients related to the rotation. In this 
work, as well as in Grids I, we chose Dii as in Zahn (1992) and 
Dshem- as in Maeder (1997). This choice relies mainly on it is 
allowing the following two behaviours simultaneously: 

- The rotating models with an initial rotational velocity that 
we assume to be the more representative (veq, ini/vcrit - 0.4, 
based on the work of Huang & Gies 2006, see also Grids I 

for more details) have a blue loop in the mass range ~ 5 

9Mq; 

- The mixing at the edge of the convective core during the MS 
is strong enough, leading to a bigger core than in the non- 
rotating case. This allows rotation to increase the size of the 
convective core in agreement with recent asteroseismic ob- 
servations by Neiner et al. (2012), who support the view that 
rotation enlarges the convective core. 

The efficiency of the mixing (the /energ factor in eq. (4) in 
Grids I) was calibrated to reproduce typical chemical enrich- 
ments at the surface of solar-metallicity MS B-type stars for our 
assumed typical initial rotation velocity (the rotating models in 
Grids I, see their Fig. 11). To be consistent with the previous 
work, we kept the same value for all the models presented here. 

3.3. Eiectronic data 

All the electronic tables of these models are available at 
http : //obswww. unige . ch/Recherche/evol/-Database-, 
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Fig. 2. MS lifetime enhancement as a function of the initial mass for all the models (same colour code as in Fig. 1). Models at Z = 0.014 (left), 
Z = 0.006 (centre), and Z = 0.002 (right). 



as well as at the CDS. We have developed an interac- 
tive web application allowing for the downloading of 
a stellar model with given Z, M, and Q/Qcrit interpo- 
lated between the computed models presented here. It 
also offers isochrones computation at the desired age 
or age range. This application is available at the url 
http : //obswww . unige . ch/Recherche/evoldb/index/. 
The database contains the models presented in this paper, as well 
as the models from Ekstrom et al. (2012) and Mowlavi et al. 
(2012). In the near future, it will be extended with additional 
metallicities. 



4. Results 

The overall characteristics at the end of the main burning stages 
for all our models are presented in Tables 2 to 4. General 
Hertzsprung-Russell diagrams (HRD) with a colour scale indi- 
cating a> are presented in Figs. 10 to 12. The shaded area repre- 
sents the approximative position of the Cepheid instability strip, 
according to Tammann et al. (2003). 

4.1. Mixing efficiency 

In Fig. 3, we show, for two masses and two metallicities, the 
profile of the diffusion coefficients responsible for the transport 
of the chemical species: the effective diffusion coefficient Dgff 
accounting for the composed effects of the meridional circula- 
tion and of the strong horizontal turbulence (see Chaboyer & 
Zahn 1992), and the shear-mixing diffusion coefficient Dshear, at 
roughly the middle of the MS phase (Xh, cen - 0.3). We see that 
except near the edge of the convective core, where Deff dom- 
inates, the shear mixing is responsible for the transport of the 
chemical species in most of the radiative envelope. We see a 
clear trend to comparing the two different masses for the same 
metallicity: both Dgff and Dsheai are higher for a higher initial 
mass. In contrast, the coefficients do not show any noticeable 
difference with respect to the metallicity. However, due to the 
lower opacities at lower metallicity, low-Z stars are more com- 
pact, leading to a shorter characteristic timescale for the diffu- 
sion (Tdiff ~ 75"), and thus a more efficient enrichment of the 



surface (see the discussion in Sect. 4.5). This behaviour remains 
true for all the models at the three considered metallicities. 

4.2. HRD, lifetimes 

As described by Meynet & Maeder (2000) and Ekstrom et al. 
(2008), the position on the ZAMS is governed by the centrifu- 
gal force: rotating models behave like lower mass ones, with a 
shift in the tracks toward lower L and Tetf. When the evolution 
proceeds, rotational mixing brings fresh hydrogen into the core, 
slowing down its decrease in mass. Also, newly produced he- 
lium is brought into the radiative zone, so the model evolves on 
a bluer and more luminous track. 

As shown in Fig. 1, the luminosity attained at the end of 
the MS, and thus during the crossing of the Hertzsprung-Russell 
gap, increases with the rotation rate up to Wini < 0.5. Beyond, the 
luminosity decreases again and the crossing occurs at a luminos- 
ity that is lower than the one of the tLtj,,; = 0.3 model, except at 
Z = 0.002, where the luminosity at the end of the MS increases 
again above Wini ~ 0.8. This behaviour is directly related to the 
size of the convective core, which governs the luminosity at the 
end of the MS. During the MS, the evolution of the size of the 
convective core relies on two counteracting physical processes 
linked to rotation. First, rotation generates an additional support 
against gravity due to the centrifugal force. These effects tend to 
decrease the size of the core and its luminosity. Second, the ro- 
tational mixing at the edge of the core progressively brings fresh 
material into the core, increasing its mass, hence its luminosity. 
Both effects are responsible for the non-monotonic behaviour of 
the luminosity at the HRD-crossing as a function of the initial 
rotation rate. 

For the models going through a blue loop during core He- 
burning (passing through a classical Cepheid phase), we note 
that rotation increases the luminosity at which the loop occurs. 
The effect of rotation on the blue loops becomes spectacular at 
low metallicity, with the disappearance of the typical loop pat- 
tern for the most rapid rotators. We discuss this point in more 
detail in Sect. 5.1. 

Rotational mixing increases the MS lifetimes, as shown in 
Fig. 2. At solar metallicity, for the mass range considered, the 
increase does not much vary as a function of the initial mass. It 
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Fig. 3. Diffusion coefficients profiles for a 2 (red curves) and a 12 Mq model (blue curves) at Z = Zq (left panel) and Z = Zsmc (right panel). The 
effective diffusion coefficient D^ff (accounting for the effects of the meridional circulation on the chemical species) is plotted in dashed lines, and 
the shear mixing diffusion coefficient in solid lines. 
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Fig. 4. MS duration enhancement as a function of the initial rotation rate. Left: 1.7 Mq models. Right: 15 Mq models. 



amounts to about 15 - 25% for Wini between 0.30 and 0.50. For 
the most rapid rotators considered here, it may reach values as 
high as 30 - 35%. 

When the metallicity decreases, the variation with the initial 
mass becomes more marked at the low and high mass ranges 
considered. As a numerical example, at Z = 0.002, the maxi- 
mum increase reaches values as high as 62% for the 1.7 Mq and 
as high as 41% for the 15 Mq. 

The high sensitivity at 1 .7 Mq comes from the central tem- 
perature at the beginning of the core H-burning phase being 
around 20 MK, i.e. in the range of temperature at which CNO- 



burning begins to dominate the pp-chain for the energy produc- 
tion. When rotation is fast enough, the hydrostatic effects of rota- 
tion decrease the central temperature, and therefore the pp-chain 
becomes the dominant energy production channel. Typically, 
at 200 Myr, the three more rapidly rotating models have < 
21.4 MK, while all the others have > 21.7MK. As a result, 
the convective core is smaller during the first part of the MS 
phase (see Fig. 5, left)*. It is amazing to note that such a slight 
difference in is enough to pass from pp-chain to CNO-cycle 

The decrease shown in Fig. 5 (left) is of the same order of magni- 
tude as the one obtained by Mowlavi et al. (2012) when the initial mass 
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dominated energy production. When evolution proceeds, and the 
central temperature increases, the convective core grows and 
reaches values even higher than the values obtained at the be- 
ginning of the core H-buming phase in the non-rotating model. 

At low metallicities and masses above 5-7 Mq, we obtain 
large increases of the MS lifetimes for Wini > 0.70. Actually, 
some small convective zones develop in the region above the 
core, where there is a variable chemical composition. (This phe- 
nomenon is related to the phenomenon of "semi-convection".) 
These convective zones mix the region above the core (and, in 
some cases, even merge with the core), bringing fresh hydrogen 
close to it, favouring a refuelling of the core, and in turn increas- 
ing the lifetime on the MS. The details in the development of 
these intermediate convective zones are very sensitive to the pre- 
cise prescription used for rotation, however the general features 
seem robust and occur independently of the physical ingredients. 

In Fig. 4, we compare the degree of dispersion of the MS 
lifetimes due to change of metallicity and of rotation for two 
initial masses. We see that for the 15 Mq (Fig. 4, right), up to 
ojjni = 0.60, the dispersion due to rotation is much greater than 
the one due to metallicity. This is due to the change in metallicity 
for the hot massive stars that have only moderate effects on the 
opacities^, hence on the luminosity of the stars. In that case, the 
mixing is the dominant effect changing the MS duration. At high 
velocities, however, the change brought by varying the metallic- 
ity are more important, but this is again a consequence of rota- 
tional mixing whose consequences are not the same depending 
on the initial metallicity (see above). 

In the low mass-range (1.7 Mq, Fig. 4, left), we somehow 
have an inverse situation, in the sense that, in general, the effects 
of changing the metallicity have a greater impact on the MS life- 
time than changing the initial velocity. This reflects the fact that, 
on the one hand, rotational mixing becomes less efficient in low- 
mass stars (see Fig. 8), and on the other hand, the impact of a 
change in Z on the opacity is much greater than in the higher 
mass range. 

4.3. Core mass vs rotation and metallicity 

As mentioned above, the additional support against gravity pro- 
duced by the rotation allows the central temperatures to be lower 
for increasing initial velocities. This implies that the size of the 
convective core on the ZAMS decreases as a function of the ini- 
tial rotation velocity. However, as evolution proceeds, the more 
efficient mixing in the more rapidly rotating stars increases the 
size of the convective core, which becomes bigger than the core 
of more slowly rotating stars. At the end of the MS, the faster 
the stars initial rotation, the larger the convective core. This be- 
haviour is valid at the three considered metallicities. We recall, 
however, the peculiar behaviour of the convective core of the 
1.7 Mq at different rotational velocities, owing to the conflicts 
between the pp-chains and CNO cycle (see above and Fig. 5, 
left panel). 

At low metalhcity, for the upper mass range of our sam- 
ple (above ~ 7 Mq), and for the fastest rotational velocities 
(cij > 0.8), the merger of the core with the convective zone that 
develops above it modifies the behavioiu' at the end of the MS, 
refuelhng it in fresh hydrogen and prolonging the MS duration 



decreases from 1.7 down to 1.6 Mq, precisely for the same effect as on 
the energy production channel. 

' The opacity in the temperature range of the massive stars is dom- 
inated by electron scattering, which is only marginally affected by a 
change in the chemical composition. 



(see Fig 5, right panel). The core mass will have an effect on the 
final yields, as is discussed in Sect. 5.2. 

4.4. Evolution of the surface velocity 

Figure 6 shows the evolution of the surface velocities during the 
MS phase for the 1.7 and 15 Mq stellar models. Close to the 
ZAMS, we see an abrupt decrease in Vgq. It corresponds to the 
time taken by the model to reach a quasi-equilibrium state (the 
initial flat Q-profile is not at the equihbrium once the meridional 
circulation acts). We also see that the equatorial velocities are 
initially higher at lower metallicity for the models starting at the 
same oj. This is due to the larger compactness of the star at low 
Z: since the initial CNO content is lower, the star must compen- 
sate for this lack of catalyst elements during the core hydrogen- 
burning phase by adopting a more compact structure that allows 
a higher central temperature. With Veq, ini = ^imf, and looking 
at eq. 1, we find that for a given co and a given mass, the initial 
equatorial velocity is Veq, ini ~ r~^, and thus, is higher for a more 
compact star. 

The surface angular velocity varies as a function of the time 
because of three main physical processes: 

- the local conservation of the angular momentum, which 
modifies the angular velocity when the star contracts or ex- 
pands; 

- internal transport mechanisms (see Zahn 1992; Maeder 
& Zahn 1998), which redistribute the angular momentum 
throughout the star; 

- stellar winds, which remove angular momentum from the 
stellar surface. 

During the MS, our models develop an external meridional cir- 
culation cell (Gratton-Opik cell, see e.g. Maeder 2009), which 
carries angular momentum from the inner part of the star to the 
surface, tending to accelerate it. The efficiency of the merid- 
ional circulation is greater for higher masses, for higher initial 
rotational rates and for higher metallicity. In contrast, the stel- 
lar mass loss due to radiative winds (Castor et al. 1975) tends to 
brake the stellar surface. As for the meridional circulation, the 
strength of the stellar winds is greater for higher mass stars, for 
higher rotational rates (Maeder & Meynet 2000), and for higher 
metallicities. There are thus two counteracting effects governing 
the evolution of the stellar surface velocity. For a given initial 
mass, rotation rate, and metallicity, their relative efficiency will 
govern the increase or decrease in the surface velocity during 
stellar fife. We see from Fig. 6 that in the mass range studied 
here, the equatorial velocity remains roughly constant during the 
MS, however, according to Eq. 3, Vcrit itself decreases during the 
MS, since /?eq is increasing. As a result, the ratio Veq/vcrit in- 
creases (see Ekstrom et al. 2008). 

From Fig. 7, one sees that the contrast between the rotation 
rates of the core and of the envelope is stronger in slowly ro- 
tating stars than in faster rotating ones. This illustrates that the 
combined effects of meridional circulation and shears, which 
tend to flatten the Q-profiles, are stronger in faster rotating stars. 
Similar qualitative behaviours are obtained for the whole mass 
and metalhcity range considered here. The model with cjim = 
0.50 (which is close to the peak of the rotation rate distribution 
according to Huang et al. 2010) shows a ratio t^cen/^^surf - 3. At 
the moment, the ratio of Ocen/^^surf has been estimated through 
asteroseismology in three slowly rotating B-type stars (Aerts 
2008). The results range between 1 and 6. The present results 
overlap the observed range well, although the small amount of 
data does not allow a systematic and detailed comparison. 
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Fig. 5. Core mass evolution. Left: 1.7 Mq models at Z = 0.002 (same colour code as in Fig. I). Right: 15 Mq models at the three Z, and fl/ficrit = 
and 0.95. 
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Fig. 6. Equatorial velocity evolution (same colour code as in Fig. 1) at Z = 0.014 (solid lines) and Z = 0.002 (dotted lines). Left: 1.7 Mq models. 
Right: 15 Mq models. 



4.5. Surface abundances 

The variations in the surface abundances at the end of the MS 
phase are shown in Fig. 8. We see, as already obtained in pre- 
vious works (Maeder & Meynet 2001), that for a given value of 
Wini*", the enrichment is stronger at lower metallicities and for 
higher masses. This is true not only for the massive stars, but 
also for intermediate-mass stars. 



* We focus our discussion on comparing models with the same Wjn;. 
The trends can differ slightly if we considered for example the same 
Vcq, ini- The detailed trends can be obtained by the mean of the tables 
provided in this paper, or directly from the available electronic data. 



In Fig. 8, the shadowed area corresponds to a variation 
of 0.2 dex with respect to the initial N/H ratio (corresponding 
roughly to the typical error bars in the measurements of abun- 
dances. For example, typical error bars in the VLT-FLAMES 
survey of massive stars for individual stars are estimated by the 
authors to be between 0.1 and 0.3 dex (Hunter et al. 2007). 
According to Nieva & Przybilla (2010), the typical error bars 
are systematically underestimated, and should rather be of the 
order of 0.3 dex. The models evolving out of this area will there- 
fore have a strong enough enrichment to be observable. The 
dashed lines correspond to the mid-MS (when the central hydro- 
gen mass fraction is half of its initial value), and the continuous 
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Fig. 7. ncentrc/^^surfacc Tatio for the 15 Mq models et Z = 0.014 (same 
colour code as in Fig. 1). 

lines to the end of the MS. At solar metallicity, the mixing is 
efficient enough in massive stars (M > 9 Mq) to be observation- 
ally detectable, even for moderately rotating (oj > 0.5) stars, and 
already at the middle of the MS. For the lower masses, this is no 
longer true, and only stars with an initial mass above 4 Mq are 
mixed enough to produce an observable enrichment at the mid- 
dle of the MS, even for the fastest rotators. At the end of the MS, 
the enrichment is observable only for the most rapidly rotating 
stars for our 1.7 Mq models (oj > 0.8). 

At Z = 0.002, we see that the more efficient mixing makes 
the surface enrichment much more easily detectable, even for 
our lowest mass models. Indeed, all the stars with initial velocity 
Wini > 0.60 (yellow tracks) show sufficient enrichment at mid- 
MS, and ojini > 0.50 for the end of the MS. 

5. Advanced phases 

5.1. Blue loops 

The physical mechanisms leading to the expansion or contrac- 
tion of the envelope during the stellar life are still poorly under- 
stood, and are being debated (Renzini et al. 1992; Stancliffe et al. 
2009). According to Lauterborn et al. (1971), the occurrence of 
a blue loop during core He-burning depends on the relation be- 
tween the gravitational potential of the core and a critical poten- 
tial that is mass dependent, as well as on the internal hydrogen 
profile. Any process able to modify either the total mass, the 
mass or radius of the core, or the hydrogen abundance profile 
will affect the occurrence of a blue loop. 

Rotation affects all these quantities in various ways, depend- 
ing on the rotation rate, as can be seen in Fig. 9. Metallicity also 
plays a role, since rotational mixing occurs differently at differ- 
ent metallicities. A general feature is that rotation increases the 
time spent on the loop. An exception to that are the models at Zq 
rotating with a>i„i - 0.8 and 0.9, where the loop is suppressed. 

At non-solar metallicity, the difference in luminosity be- 
tween the leftwards excursion at the start of the loop and the 
redwards return movement is widened by rotation. The widen- 



ing occurs toward both lower starting luminosity and higher final 
luminosity. In the most extreme cases, the lowering of the start- 
ing luminosity might reach the luminosity of the crossing of the 
Hertzsprung gap. In that case, the model "jumps" directly onto 
the top of the loop, avoiding the first RSG phase. For those mod- 
els, most of the core He-burning occurs thus in the blue part of 
the HRD (log(reff [K]) > 4.0). Such behaviour will reduce the 
number of RSG at low metallicity. 

We will dedicate a future paper to more detailed study of the 
implication of this effect on the Cepheid population predicted by 
our models. 

5.2. Mass ejections of various elements 

In Table 1, we present the ejected masses for H, He, C, N, O, 
and remaining metals, as well as the CO-core mass and the rem- 
nant mass (computed as in Hirschi et al. 2005, and Georgy et al. 
2009). We see that rotation increases the mass of the CO core 
significantly, and this in a larger way at lower metallicity, which 
matches the result of a higher mixing in lower Z. For the SMC, 
the increase amounts to 43%, while for the LMC and Galaxy it 
is 27% and 25%, respectively, for the models with Wi„i - 0.95. 
This table is qualitatively representative of the trends obtained 
for the higher part of our mass range. We do not discuss here 
the ejected masses for the models of our sample that become 
AGB stars, since we have not followed these models through 
this phase. 

The general trends concerning the ejected elements in rotat- 
ing models are the following. 

- Hydrogen is mostly depleted (-60% in the SMC and LMC 
models, and 73% in the Galactic models); this reflects the 
fact that more mass is processed by nuclear burning in rotat- 
ing stars. 

- Carbon and oxygen tend to be produced more. Comparing 
the models with Wini = and Wini - 0.95, the increase in the 
ejected mass of carbon amounts to a factor of 2.15 for the 
SMC models, 1.52 for the LMC, and 1.71 for the Galactic 
models. The increase in the ejected mass of oxygen amounts 
to a factor of 3.36 for the SMC, 2.52 for the LMC, and 2.40 
for the Galaxy. In the Galactic case, the higher increase in C 
and O occurs for models with Winj - 0.50 and 0.60, respec- 
tively, with a maximal increase by a factor of 2.06 and 3.11, 
respectively. 

Helium, nitrogen and the remaining metals do not show any def- 
inite trend. One also notes that even at high rotation, we do not 
observe any production of primary nitrogen, in these models. 

6. Conclusions 

In this paper, we have presented an extended database of rotat- 
ing stellar models at three different metallicities, for nine dif- 
ferent initial rotation parameters and ten different masses, in 
the line of the previous large grid of stellar models (Ekstrom 
et al. 2012). The computations account accurate following of the 
angular-momentum content, and for the stellar-wind anisotropy. 
Moreover, they allow to give an estimate of the amount of mass 
that the star should lose mechanically in an equatorial disc when 
rotating at the critical velocity. This database will be particularly 
useful for constructing synthetic populations of stars, accounting 
for mass, rotation, and metallicity distributions. 

In this first paper, we presented some general results on the 
behaviour of the models in the HRD, their lifetime, surface ve- 
locities, chemical enrichment, and advanced phases behaviour. 
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Fig. 8. N/H ratio at the end of the MS as a function of the initial mass for all the models. Models at Z = 0.014 {left), Z = 0.006 (centre), and 
Z = 0.002 (right). 
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Fig. 9. Teff as a function of the helium mass fraction in the core. Models at Z = 0.014 (left), Z = 0.006 (centre), and Z = 0.002 (right). 



In the second paper of this series (Granada et al. 2012), extensive 
comparison between the results presented here and observation 
concerning Be stars were done, particularly concerning the mean 
mass-loss rates during the critical-rotation phase, disc mass, and 
lifetimes. 
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Table 1. CO-core mass, remnant mass, and ejected masses in hydrogen, helium, carbon, nitrogen, oxygen, and remaining heavy elements for the 
15 Mq models at the three metallicities and all the rotation rates. All masses are expressed in M© units. 





Mco 






^He 


Mc 






M 

other heavies 


Z = 0.002 0.00 


2.6961 


1.4265 


6.8329 


5.0423 


0.3638 


0.0056 


0.6237 


0.2129 


0.10 


3.0113 


1.4915 


6.5290 


4.8647 


0.3795 


0.0058 


0.9746 


0.1920 


0.30 


3.1037 


1.5097 


6.0193 


5.1801 


0.4782 


0.0066 


1.1124 


0.1496 


0.50 


3.0275 


1.4947 


5.5969 


5.2227 


0.5249 


0.0071 


1.1904 


0.1022 


0.60 


2.8358 


1.4556 


5.8074 


5.5233 


0.5636 


0.0084 


1.0755 


0.0744 


0.70 


3.3706 


1.5619 


5.1903 


4.8331 


0.5772 


0.0071 


1.5936 


0.1491 


0.80 


3.3276 


1.5535 


4.9495 


5.2080 


0.6117 


0.0079 


1.5929 


0.1397 


0.90 


3.7427 


1.6335 


4.3507 


4.8306 


0.6983 


0.0069 


1.9760 


0.0207 


0.95 


3.8653 


1.6568 


4.1340 


4.7515 


0.7848 


0.0069 


2.0926 


0.0442 



Z = 0.006 


0.00 


2.6607 


1.4191 


6.2795 


4.8178 


0.4218 


0.0167 


0.6678 


0.1142 




0.10 


2.8571 


1.4600 


5.7171 


4.5974 


0.3637 


0.0162 


0.8697 


0.2075 




0.30 


3.1193 


1.5128 


5.2323 


4.7427 


0.4892 


0.0167 


1.1535 


0.1483 




0.50 


3.0771 


1.5045 


5.2111 


5.0265 


0.5610 


0.0185 


1.1342 


0.1446 




0.60 


2.9607 


1.4815 


5.1949 


5.0690 


0.5576 


0.0192 


1.1200 


0.1409 




0.70 


3.1077 


1.5105 


5.2193 


4.7928 


0.6010 


0.0203 


1.4027 


0.1191 




0.80 


3.1568 


1.5202 


4.9762 


4.7993 


0.5864 


0.0204 


1.5126 


0.1071 




0.90 


3.4196 


1.5714 


3.9929 


4.4800 


0.6119 


0.0180 


1.6830 


0.1459 




0.95 


3.4000 


1.5676 


3.8333 


4.3666 


0.6419 


0.0175 


1.6814 


0.1237 


Z = 0.014 


0.00 


2.4065 


1.3655 


5.8070 


4.8193 


0.2956 


0.0400 


0.4983 


0.2577 




0.10 


2.7714 


1.4422 


4.8298 


4.2007 


0.3913 


0.0337 


0.8906 


0.1801 




0.30 


3.3153 


1.5511 


3.4444 


3.5356 


0.5674 


0.0272 


1.4047 


0.1666 




0.50 


3.1867 


1.5260 


3.4473 


3.7215 


0.6409 


0.0287 


1.3235 


0.1908 




0.60 


3.2256 


1.5337 


3.4455 


3.8259 


0.5787 


0.0309 


1.4612 


0.1081 




0.70 


3.2558 


1.5395 


3.4518 


3.9198 


0.5009 


0.0334 


1.5509 


0.0936 




0.80 


2.9022 


1.4694 


4.1405 


4.7817 


0.5307 


0.0392 


1.0723 


0.1649 




0.90 


2.9977 


1.4888 


4.0481 


4.3689 


0.6091 


0.0382 


1.2925 


0.1168 




0.95 


3.0261 


1.4944 


4.2465 


4.6714 


0.5045 


0.0420 


1.1951 


0.1783 



Notes. The remnant mass is computed following the same method as in Hirschi et al. (2005) and Georgy et al. (2009). 
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Table 2. Main parameters of A-B stars at Z = 0.014 = Z© 



Mo 



End of H-burning 



N/C N/0 
lass fract. 



End of He-burning 



End of C-buming 





M 








I'surf 


N/C 


N/0 


kyr 


Mo 


day 








mass fract. 




4.721 


13.25 










0.3390 


3.0660 


0.7360 


2.868 


12.14 


3.780 


10* 


3.192 


10^ 


0.3480 


3.0057 


0.7913 


1.516 


10.89 


1.429 


10^ 


2.535 


10^ 


0.3769 


4.3558 


0.9605 


2.121 


11.07 


1.106 


10'" 


2.685 


10" 


0.3954 


6.5423 


1.1182 


1.540 


11.17 


7.067 


10^ 


1.974 


10" 


0.4161 


8.8192 


1.2717 


1.216 


11.29 


3.946 


10' 


1.142 


lO" 


0.4264 


11.7377 


1.3780 


2.687 


12.36 


1,263 


lO' 


4.776 


lO*' 


0.4207 


13.9606 


1.3818 


2.479 


12.13 


1,445 


105 


6,067 


10« 


0.4180 


16.9923 


1.4089 


2.543 


12.51 


9.546 


Itf 


4.241 


10^ 


0.4208 


20.0816 


1.4607 


7.858 


11.31 










0.3072 


1 .9788 


0.5280 


6.859 


11.15 


4.957 


10' 


6.823 


10^ 


0.3082 


2.4007 


0.6104 


6.705 


10.57 


2.214 


10' 


6.239 


10^ 


0.3287 


3.2817 


0.7213 


6.134 


10.31 


1.327 


10' 


5.280 


lO'' 


0.3502 


5.3748 


0.8895 


3.803 


10.22 


1.292 


10^ 


5.572 


10* 


0.3611 


6.9967 


0.9729 


6.097 


10.76 


5,262 


10* 


3,531 




0-3722 


9.6176 


1.0656 


0.063 


10.64 


5,642 


10* 


5,408 




0.3712 


12.7958 


1.1021 


0.072 


10.69 


5,110 


10* 


5,208 


10' 


0.3746 


17.8632 


1,1482 


0.066 


10.73 


4,571 


10* 


5,457 


10' 


0.3780 


21.0351 


1,1874 



15,00 
15,00 

15.00 
15.00 
15.00 
15.00 
15.00 
15.00 
15.00 
12.00 
12.00 
12.00 
12.00 
12.00 
12.00 
12.00 
12.00 
12.00 
9.00 
9.00 
9.00 
9.00 
9.00 
9.00 
9.00 
9.00 
9.00 
7.00 
7.00 
7.00 
7,00 
7,00 
7.00 
7.00 
7.00 
7.00 
5.00 
5.00 
5,00 
5,00 
5.00 
5.00 
5.00 
5.00 
5.00 
4.00 
4.00 
4.00 
4.00 
4.00 
4.00 
4.00 
4.00 
4.00 
3.00 
3.00 
3.00 
3.00 
3.00 
3.00 
3.00 
3,00 
3,00 
2.50 
2.50 
2.50 
2.50 
2.50 
2.50 
2.50 
2.50 
2.50 
2.00 
2.00 
2.00 
2.00 
2.00 
2.00 
2.00 
2.00 
2.00 



aoo 

0.10 
0.30 
0.50 
0.60 
0.70 
0.80 
0.90 
0.95 
0.00 
0.10 
0.30 
0.50 
0.60 
0.70 
0.80 
0.90 
0.95 
0.00 
0.10 
0.30 
0.50 
0.60 
0.70 



0.00 
0.10 
0.30 
0.50 
0.60 
0.70 
0.80 
0.90 
0.95 
0.00 
0.10 
0.30 
0.50 
0.60 
0.70 
0.80 
0.90 
0.95 
0.00 
0.10 
0.30 
0.50 
0.60 
0.70 



0.00 
0.10 
0.30 
0.50 
0.60 
0.70 
0.80 
0.90 

_a95_ 

0.00 
0.10 
0.30 
0.50 
0.60 
0.70 
0.80 
0.90 
0.95 
0.00 
0.10 
0.30 
0.50 
0.60 
0.70 
0.80 
0.90 
0.95 
0.00 
0.10 
0.30 
0.50 
0.60 
0.70 
0.80 
0.90 
0.95 



11.015 

II. 661 

12.842 
13.399 
13.633 
13.836 
13.983 
14.123 
14.189 
15.331 
16.182 
17.862 
18.540 
18.762 
18.902 
19.000 
19.126 
19.302 
26.262 
27.541 
30.354 
31.464 
31.685 
31.918 
31.971 
32.264 
32.482 
41.722 
44.65 1 
49.304 
51,151 
51.688 
51.999 
52.238 
52.466 
52.900 
88.195 
94.322 
104.456 
109.085 
110.343 

III. 232 
112.264 
112.853 
113.858 
152.088 
160.880 
178.575 
187.477 
190.190 
192.706 
194.065 
195.595 
198.351 
320.600 
338.503 
378.862 
400.944 
408.242 
413.122 
417.548 
420.785 
422.412 
537.964 
555.638 
622.923 
665.370 
677.421 
687.051 
695.645 
702.516 
704.202 
1008.888 
1050.551 
1177.056 
1267.754 
1293.495 
1316.637 
1332.516 
1337.360 
1340.230 
1633.971 
1705.879 
1903.093 
2068.535 
2108.224 
2148.808 
2166,672 
2200.253 
2181.676 



14.81 
14.79 

14.73 



8.99 
8.99 
8.98 



7.00 
7.00 
7.00 
7.00 
7.00 
7.00 
6.99 
5.00 
5.00 
5.00 
5.00 
5.00 
5.00 
5.00 
5.00 



3.00 
3.00 
3.00 
3.00 
3.00 
3.00 
3.00 
3.00 
3.00 
2.50 
2.50 
2.50 
2.50 
2.50 
2.50 
2.50 
2.50 
2.50 
2.00 
2.00 
2.00 
2.00 
2.00 
2.00 
2.00 



1.70 
1.70 
1.70 
1.70 



0.000 

0.043 
O.I 82 
0.329 
0.409 
0.480 
0.550 
0.614 
0.639 
0.000 
0.054 
0.216 
0.397 
0.493 
0.594 
0.701 
0.749 
0.748 
0.000 
0.060 
0.230 
0.427 
0.527 
0.632 
0.758 
0.887 
0.888 
0.000 
0.060 
0.226 
0.420 
0.527 
0.641 
0.788 
0.905 
0.892 
0.000 
0.060 
0.218 
0.409 
0.516 
0.629 
0.783 
0.904 
0.866 
0.000 
0.058 
0.210 
0.399 
0.506 
0.627 
0.795 
0.849 
0.849 
0.000 
0.057 
0.205 
0.395 
0.497 
0.614 
0.788 
0.997 
0.748 
0.000 
0.058 
0.200 
0.389 
0.496 
0.616 
0.784 
0.827 
0.800 
0.000 
0.060 
0.203 
0.392 
0.502 
0.629 
0.799 
0.826 
0.775 
O.(K)0 
0.062 
0.202 
0.391 
0.498 
0.625 
0.783 
0.821 



0,2660 
0.2660 

0.2704 
0.2835 
0.2922 
0.3008 
0.3087 
0,3164 
0.3199 
0.2660 
0.2660 
0.2674 
0.2755 
0,2813 
0,2879 
0,2944 
0.3026 
0.3071 
0.2660 
0.2660 
0.2664 
0.2709 
0.2748 
0.2799 
0.2847 
0.2913 
0.2962 
0.2660 
0.2660 
0.2661 
0,2685 
0,2711 
0.2743 
0.2774 
0.2826 
0.2871 
0.2660 
0,2660 
0,2660 
0.2670 
0.2684 
0.2703 
0.2715 
0.2733 
0.2771 
0.2660 
0.2660 
0.2660 
0.2666 
0.2674 
0.2687 
0.2696 
0.2712 
0.2743 
0.2660 
0.2660 
0.2660 
0.2664 
0.2670 
0.2677 
0.2680 
0,2697 
0,2719 
0.2660 
0.2660 
0.2660 
0.2663 
0.2668 
0,2674 
0,2676 
0.2691 
0.2711 
0.2660 
0.2660 
0.2660 
0.2664 
0.2669 
0.2675 
0.2677 
0.2696 
0.2704 
0.2660 
0.2660 
0.2660 
0.2664 
0.2670 
0,2676 
0.2677 
0,2694 
0,2711 



0.2885 
0.2931 

0.7316 
1.7907 
2.6259 
3.7090 
5.0934 
6.9790 
8.1555 
0.2885 
0.2886 
0.5172 
1.4363 
2.1637 
3.2498 
4.7996 
7.3874 
9.0862 
0.2885 
0.2885 
0.3841 
1.0243 
1.5832 
2.4400 
3.6168 
5.4922 
7.0486 
0.2885 
0.2885 
0.3302 
0.7166 
1.1053 
1.6250 
2.2416 
3.3200 
4.3773 
0.2885 
0.2885 
0.3026 
0.4995 
0.7165 
1.0124 
1.2602 
1.5533 
2.1374 
0.2885 
0.2885 
0.2951 
0.4202 
0.5617 
0.7616 
0.9224 
1.1427 
1.5901 
0.2885 
0.2885 
0.2913 



0.5880 
0.6515 
0.8674 
1.1641 
0.2885 
0.2885 
0.2900 
0.3491 
0.4262 
0.5150 
0.5592 
0.7478 
1.0004 
0.2885 
0.2885 
0.2892 
0.3352 
0.3978 
0.4682 
0.4900 
0.6955 



0.1152 
0.1166 

0.2442 
04568 
0.5704 
0.6758 
0.7694 
0.8560 
0-8963 
0.1152 
0.1152 
0.1846 
0-3818 
0.4829 
0.5875 
0.6845 
0.7816 
0.8352 
0.1152 
0.1152 
0.1454 
0.3023 
0.3960 
0.4972 
0-5868 
0.6907 
0.7520 
0.1152 
0.1152 
0-1285 
0.2341 
0.3149 
0-3975 
0.4701 
0.5621 
0.6282 
0-1152 
0.1152 
0.1197 
0-1781 
0.2329 
0.2950 
0.3371 
0-3806 
0.4527 



0.1172 
0.1554 
0.1941 
0-2421 
0.2750 
0.3159 
0.3846 
0-1152 
0-1152 
0.1161 



0.1993 
0.2140 
0.2608 
0.3145 
0-1152 
0.1152 
0.1156 
0.1338 
0-1559 
0.1795 
0.1903 
0-2333 
0.2816 
0.1152 
0.1152 
0.1154 
0-1294 
0-1474 
0.1663 
0.1718 
0.2187 



0.7895 0.2373 

0.2885 0.1152 

0.2885 0.1152 

0.2887 0.1153 

0.3214 0.1252 

0.3732 0-1402 

0.4319 0.1562 

0.4483 0.1605 

0.5730 0.1907 

0.7294 0.2228 



1.441 

2.124 
2.150 
1.812 
1.892 
2,125 
2,006 
2,178 
2,400 
2,357 
3,492 
3,502 
3.296 



3,450 
3,610 
4,446 
4,356 
4,610 
6,919 
7,138 
6,663 
6,795 
7,174 
6,990 
7,431 
8,239 
8,465 
19,380 
18.712 
17.918 
16.800 
17.669 
17.447 
18.920 
20.244 
21.160 
37.750 
42.425 
35.069 
34.852 
33.704 
33.980 
35.455 
37.172 
37.370 



111.205 
97.083 
98.701 
101.601 
102.865 
109.725 
105.941 
236-764 
253.558 
205.343 
185.800 
191-139 
182.884 
193.384 
193.248 
194.550 



1141 
12.48 
12.26 
12.63 



2.833 
1.424 
7.856 
4.914 
2.126 
7.165 
9.525 



10.65 
10.40 
10.30 
10.84 
10.71 
10.76 
10.80 



8.78 
8.78 



6.92 
6.92 
6.90 



6.89 
6.90 
6.90 
6.89 
4.96 
4,96 
4,95 
4,95 
4,95 
4,95 
4,95 
4,95 
4,94 



3,97 
3,97 
3,97 
3,97 
3,97 
3,96 



2,98 
2,98 
2,98 



2,49 
2,49 
2,49 
2,49 
2,48 
2,49 
2,48 
2.48 
2.48 



1.665 
1.205 
8.402 
7.799 
3.031 
3.647 
3.236 
3.044 



2.852 
1.501 
l.(H)l 
8.019 
6.505 
6.841 
6.055 
9.626 

1.679 
6.959 
4.772 
4.221 
3.169 
2.645 
2.440 
2.508 



4.324 
3.897 
3.954 



1.577 
1.471 
1.478 



2.579 
2.842 
2.153 
1.583 
7.387 
2.802 
4.267 
2.815 



3.027 
1.370 



3.767 
3.524 
4.341 



2..? 

2.238 

2.193 

2.268 

2.037 

1.956 

2.059 

2.187 



1.415 
9.568 
7.935 
8.018 



5.097 
4.237 
4.076 
4.007 
4.132 
4.232 

2.258 
2.201 
2.010 

2.032 
1.991 
2.061 
2.026 
2.302 



0,3054 

0,3235 

0.3564 



0,4100 
0.4022 
0,3976 
0,4024 
0,3007 
0,3053 
0,3193 
0,3272 
0,3362 
0,3518 
0,3428 
0,3508 
0,3552 
0,2824 
0,2924 
0,3169 
0,3285 
0,3330 
0,3357 
0,3334 
0,3285 
0..3296 
0.2792 
0.2880 
0,3104 
0,3190 
0,3221 
0,3248 
0,3266 
0,3285 
0,3291 
0,2792 
0.2859 
0.3071 
0.3154 
0.3190 
0.3205 
0.3225 
0.3232 
0.3256 
0.2823 
0.2873 
0..3077 



[{f 0.3187 

10^ 0.3213 

10* 0.3227 

10* 0.3245 

10* 0.3276 

0.2868 

10* 0.2906 

10* 0.3106 

10* 0.3213 

10* 0.3244 

10* 0.3261 

10* 0.3281 

10* 0.3278 

10* 0,3314 

0,2852 

10* 0,2884 

10* 0..3O80 

1(V 0,3206 

10* 0,3230 

10* 0,3249 

10* 0,3271 

10* 0,3299 

10* 0,3301 



2,0811 

2,3726 
3.7688 
5.5490 
7.6476 
10.2467 
11.9872 
14.7608 
17.4288 
1.8374 
2.3165 
3.0294 
4.4674 
5.8284 
8,1420 
10.5711 
15.2610 
18.2425 
1.6090 
2.0461 
3.0518 
4.7700 
6.1149 
8.1696 
10.8419 
13.7674 
15.9175 
1.5177 

I. 8383 
2.7399 
4.0102 
4.9466 
6.3624 
8.3759 

II. 3917 
13.0544 
1.4819 
1,6282 
2,4483 
3,3394 
4,0502 
4,8353 
5,8563 
7,1672 
8,1090 
1,5023 
1,5716 
2,2908 
3.0686 
3,5117 
4,1455 
4,9073 
5,7566 
6,4637 



3,3040 
3,7223 
4,3228 
4,8606 
5,3924 
1,4444 
1,4097 
2,0094 
2,6893 
3,0078 
3,3772 
3.8853 
4.5669 
4.7606 



0.5211 
0.6413 
0.8513 
0.9928 
1.1547 
1.2709 
1.2671 
1.2829 
1.3403 
0.4918 
0.5937 
0.6752 
0.7804 
0.8577 
0.9623 
0.9748 
1.0404 
1.0840 
0.4219 
0.5234 
0.6696 
0.7896 
0.8446 
0.8956 
0.9206 
0.9446 
0.9689 
0.3964 
0.4734 
0.6210 
0.7192 
0.7671 
0.8174 
0.8637 
0.9088 
0.9246 
0.3864 
0,4269 
0,5752 
0,6627 
0,7107 
0,7467 
0,7836 
0,8163 
0,8382 



0,5518 
0.6380 
0,6695 
0,7103 
0,7418 
0,7737 
0,7997 
0,3928 
0,3994 
0,5285 
0,6132 
0,6518 
0,6762 
0,7081 
0,7256 
0,7517 
0,3662 
0,3648 
0,4834 
0,5737 
0,5997 
0.6262 
0.6559 
0.6898 
0.6958 



5.781 ■ 10* 7.649-10' 0.3172 
3.673 ■ 10* 7.449 ■ 10' 0.3290 



3.2003 0.6840 
4.9811 0.8029 



Notes. The data presented in this Table are: initial mass (column 1), initial rotation rate (column 2), initial equatorial velocity (column 3), mean 
equatorial velocity during the MS (column 4). For each burning phase, we also give the duration of the phase (columns 5, 12, and 19), the mass 
at the end of the phase (columns 6, 13, and 20), the equatorial velocity (column 7), or the rotation period at the end of the phase (columns 14, 
and 21), the ratio of the equatorial velocity to the critical one (column 8) or the of the angular velocity to the critical one at the end of the burning 
phase (columns 15, and 22), the surface helium mass fraction (columns 9, 16, 23), the ratio N/C (columns 10, 17, and 24) and N/O at the end of 
the phase (columns 11, 18, and 25). 
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Georgy et al.: Populations of rotating stars 




4.4 4.2 4 3.8 3.6 



log (T,„ [K]) 

Fig. 11. Same as Fig. 10, but for Z = 0.006. 
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Georgy et al.: Populations of rotating stars 




4.4 4.2 4 3.8 3.6 
log (T,„ [K]) 

Fig. 12. Same as Fig. 10, but for Z = 0.002. 
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Georgy et al.: Populations of rotating stars 



Table 3. Main parameters of A-B stars at Z = 0.006 = Zlmc 



End of H-burning 



End of He-burning 



Mq 




"'km 


"MS 




M 
Mq 








N/C 
mass fract. 


N/0 




M 

Mq 




-1 




I'surf 


N/C 
mass fract. 


N/0 


15.00 
15.00 


0.00 
0.10 


0. 
48. 


0. 

33. 


11.097 
11.786 


14.87 
14.85 


0. 
22. 


0.000 
0.049 


0,2559 
0.2559 


0.2885 
0.2928 


0.1152 
0.1165 


1.338 
1.403 


14,00 
13.51 


_ 

5 All 
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5 194 


(, 
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0.3566 
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15.00 
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149. 
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14.81 


91. 
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0.3412 
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15.00 
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1 143 


10'' 


0.3490 


6.5843 


1.0305 


15.00 


0.60 


289. 


226. 


13.547 


14.80 


211. 


0.458 


0.2864 


4.2331 


0.7202 


1.307 


13.89 


3.222 


10-^ 


■ 


in'> 


0.3402 


8.2676 


1.0503 


15.00 


0.70 


369. 


270. 


13 696 


14.7') 


252. 


0.536 


0.2957 


6.4.' 50 


O.S.''J') 


1.3')9 


I. '.96 
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,, 
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12.S962 
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1 .'7 12 


1 .21 1 
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4 1)"^') 




11.4.' Id 


4') 7.^5.1 


1 .(i5'.= 


12.00 


0.00 


0. 


0. 


15.310 


11.95 


0. 


0.000 


11.255') 


0.2t485 


0.1 152 


1 .976 


10.95 


_ 








0.2667 


1.5412 


0.381 1 


12.00 


0.10 


44. 


33. 


16.277 


11.95 


24. 


0.055 


0.2559 


0.2888 


0.1153 


1.942 


10.76 


4.871 


10"' 


5 427 


10'' 


0.2756 


1 .9693 


0.5192 


12.00 


0.30 


131. 


102. 


17.783 


11.93 


95. 


0.217 


0.2584 


0.7459 


0.2469 


1.822 


10.39 


1.428 




3 447 


10'' 


0.3000 


3.4505 


0.7163 


12.00 


0.50 


242. 


178. 


18.395 


11.93 


177. 


0.400 


0.2699 


2.3935 


0.5259 


1.910 


11.42 


2.350 


10* 


8.443 


10^ 


0.3068 


5.3803 


0.8292 


12.00 


0.60 


275. 


219, 


18.572 


11.93 


221. 


0.490 


0,2779 


3.9455 


0.6631 


2.046 


11,47 


1.858 


10* 




lOS 


0.3213 


8.1782 


0.9544 


12,00 
12.00 


0.70 
0.80 


330. 
423. 


262. 
308. 


18.664 
18.758 


11.93 
11.93 


267, 
318. 


0.582 
0.675 


0J2861 
0.2939 


6.4532 
10.4655 


0.7875 
0.8913 


2.065 
2.291 


11.33 
11.33 


1.975 
1.714 


10* 
10* 


9.419 


Itf 

10' 


0.3219 
0.3331 


12,4816 
20.5658 


1.0198 
1.1150 


12.00 


0.90 


510. 


360. 


21.313 


11.90 


365. 


0.777 


0.3409 


29.6185 


1.1881 


1.872 


10.98 


2.543 


10* 


^547 




0.3892 


50.5192 


1 .4042 


12.00 


0.95 


569. 


391. 


22.249 


11.88 


378. 


0.799 


0.3638 


44.9194 


1.3331 


1.802 


1 1.07 


2.209 


10-^ 


1.317 


10'' 


0.4041 


67.7565 


1.5184 


9.00 


0.00 


0. 


0. 


25.992 


8.99 


0. 


0.000 


0.2559 


0.2885 


0.1 152 


3.363 


8.65 


_ 








0.2614 


1.51. '4 


0.3785 


9.00 


0.10 


45. 


32. 


27.333 


8.99 


25. 


0.059 


0.255'* 


0.2885 


0.1 152 


3.304 


8.73 


5.180 


10"^ 


6.417 


'"; 


0.2729 


2.1345 


0.5357 


9.00 


0.30 


123. 


98. 


29.932 


8.99 


94. 


0.225 


0.2570 


0.5544 


0.1958 


3.16K 


8.70 


1.881 


10"^ 




In' 


0.2981 


3.5358 


0.7146 


9.00 


0.50 


228. 


169. 


30.861 


8.99 


175. 


0.410 


0.264') 


l.')04l 


0.4538 


3.242 


8.74 


1.244 


10"^ 


s mi 
■ 




0.3103 


6.2920 


0.8623 


9.00 


0.60 


279. 


207. 


31.091 


8.99 


218. 


0.502 


0.2713 


3.2512 


0.5922 


3.355 


8.80 


1 .069 


10"^ 






0.3098 


8.5053 


0.9093 


9.00 


0.70 


331. 


248. 


31.183 


8.99 


265. 


0.598 


0.2782 


5.4215 


0.7187 


3.749 


8.91 


6.854 


10' 


■ 




0.3074 


11.2245 


0.9394 


9.00 


0.80 


392. 


292. 


31.292 


8.99 


318. 


0.704 


0.2852 


8.8693 


0.8277 


3.788 


8.94 


5.431 


10' 






0.3132 


17.8575 


1.0164 


9.00 


0.90 


465. 


343. 


34.282 


8.99 


393. 


0.875 


0.3093 


21.1362 


1.0484 


3.282 


8.92 


5.699 


10' 


3 973 




0.3447 


39.2807 


1.2232 


9.00 


0.95 


528. 


372. 


35.516 


8.98 


394. 


0.858 


0.32'J6 


33.6904 


1.1950 


3.504 


8.90 


6.221 


10' 


4.743 


lo^ 


0.3678 


58.1684 


1.3719 


7.00 


0.00 


0. 


0. 


41.395 


7.00 


0. 


0.000 


0.2559 


0.2885 


0.1 152 


6.275 


6.89 










0.2587 


1.3783 


0.3438 


7.00 


0.10 


43. 


31. 


43.205 


7.00 


25. 


0.060 


0.2559 


0.2885 


0.1 152 


6.197 


6.90 


2.107 


10-* 


3.040 


10^ 


0.2689 


1.8983 


0.4763 


7.00 


0.30 


129. 


93. 


47.420 


7.00 


90. 


0.223 


0.2564 


0.4296 


0.1594 


5.752 


6.87 


9.951 


lO' 


3.235 


"^s 


0.2903 


3.0605 


0.6488 


7.00 
7.00 


0.50 
0.60 


217. 
260. 


160, 
197, 


48.956 
49.298 


7.00 
7.00 


166, 
207. 


0.407 
0.499 


0.2613 
0,2661 


1.3475 
Z2769 


0.3665 
0.4948 


5.903 
6.177 


6.88 
6.90 


5.753 
4.467 


Itf 
Itf 


2.729 




0J02O 
0.3059 


5.2802 
7.2685 


0.7972 
0.8655 


7.00 
7.00 


0.70 
0.80 


320. 
370. 


235. 
278. 


49.513 
49.821 


7.00 
7.00 


252. 
304. 


0.597 
0.707 


0.2716 
0.2780 


3.6906 
5.8915 


0.6153 
0.7285 


6.791 
7.209 


6.91 
6.93 


3.653 
2.956 


Itf 
Itf 


2422 
2.222 


10' 
10' 


0.3093 
0.3113 


10.3670 
14.6980 


0.9271 

0.9800 


7.00 


0.90 


447. 


327. 


50.087 


7.00 


377. 


0.859 


0.2841 


9.1889 


0.8396 


7.361 


6.94 


2.541 




2.187 


10' 


0.3101 


19.1856 


1.0263 


7.00 


0.95 


496. 


355. 


50.622 


7.00 


399. 


0.903 


0.2898 


12.0638 


0.9159 


7.729 


6.98 


2.594 




2.244 


10 


0.3155 


23.9778 


1.0888 


5.00 


0.00 


0. 


0. 


83.169 


5.00 


0. 


0.000 


0.2559 


0.2885 


0.1152 


15.908 


4.95 










0.2582 


1.3024 


0.3256 


5.00 
5.00 


0.10 
0.30 


36. 
109. 


29. 
87. 


86.831 
95.898 


5.00 
5.00 


23, 
82. 


0.060 
0.218 


0,2559 
0.2561 


0.2885 
0.3495 


0.1152 
0.1347 


16,797 
14,949 


4.94 
4.93 


7.229 
3.656 


1^ 
Itf 


1.199 
1.426 


IC 
10' 


0.2660 
0.2854 


1.6548 
2.6598 


0.4236 
0.5938 


5.00 


0.50 


187. 


150. 


99.280 


5.00 


152. 


0.401 


0.2588 


0.8924 


0.2773 


15.334 


4.93 


2.292 


105 


1.291 




0.2958 


4.2339 


0.7274 


5.00 


0.60 


245. 


184. 


100.247 


5.00 


190. 


0.496 


0.2617 


1.4211 


0.3773 


16.006 


4.93 


1.823 


10" 


1.192 




0.2999 


5.3840 


0.7876 


5.00 


0.70 


276. 


220. 


100.700 


5.00 


232. 


0.598 


0.2653 


2.1711 


0.4758 


16.435 


4.94 


1.461 


10' 


1.141 


w 


0.3024 


6.9844 


0.8332 


5.00 


0.80 


360. 


261. 


101.328 


5.00 


283. 


0.718 


0.2695 


3.2562 


0.5746 


16.528 


4.95 


1 .073 


10"' 


9.861 




0.3070 


9.8091 


0.8985 


5.00 


0.90 


421. 


307. 


101.391 


5.00 


357. 


0.890 


0.2716 


4.2426 


0.6364 


17.932 


4.95 


8.869 


10= 


9-326 


10" 


0.3072 


12.7927 


0.9371 


5.00 


0.95 


438. 


334. 


101.879 


5.00 


358. 


0.888 


0.2728 


4.6592 


0.6626 


18.136 


4.95 


8.985 


10= 






0.3070 


14.3721 


0.9593 


4.00 


0.00 


0. 


0. 


137.775 


4.00 


0. 


0.000 


0.2559 


0.2885 


0.1152 


33.275 


3.96 






- 




0.2623 


1.4460 


0.3596 


4.00 


0.10 


38. 


28. 


143.820 


4.00 


22. 


0.059 


0.2559 


0.2885 


0.1152 


33.000 


3.96 


3.350 


10" 


6.122 


10* 


0.2682 


1.6244 


0.4137 


4.00 


0.30 


115. 


83. 


159.264 


4.00 


76. 


0.215 


0.2560 


0.3208 


0.1255 


29.021 


3.95 


1.572 


10" 


6.698 


10" 


0.2874 


2.5587 


0.5779 


4.00 


0.50 


196. 


143. 


165.639 


4.00 


141. 


0.396 


0.2576 


0.6796 


0.2261 


28.500 


3.95 


9.016 


10= 


5.691 


10" 


0.2970 


3.8659 


0.6')70 


4.00 


0.60 


240. 


176. 


167.463 


4.00 


177. 


0.494 


0.2598 


1.0704 


0.3111 


28.864 


3.95 


7.179 


10= 


5.405 


10" 


0.3001 


4.8095 


0.7495 


4.00 


0.70 


282. 


211. 


168.745 


4.00 


218. 


0.604 


0.2624 


1.5673 


0.3931 


30.519 


3.96 


6.436 


10= 


5.558 


10" 


0.3030 


5.9616 


0.7940 


4.00 
4.00 
4.00 


0.80 
0.90 
0.95 


332. 
403. 
471. 


249, 
295. 

322. 


169.522 
170.134 
170.473 


4.00 
4.00 
4.00 


268. 
308. 
302. 


0.732 
0.830 
0.810 


0.2651 
0.2650 
0.2692 


2,1962 
2.4387 

3.3039 


0.4704 
0.4900 

0.5652 


32,690 
33.696 
35.440 


3.96 
3.96 
3.96 


5.762 
4.728 
4.810 


10^ 
10^ 

10^ 


5.656 
5.373 


10* 
10* 


0.3067 
0.3089 
0.3085 


7.8889 
10.4870 
11.7103 


0.8485 
0.8932 
0.9063 


3.00 


0.00 


0. 


0. 


276.224 


3.00 


0. 


0.000 


0.2559 


0.2885 


0.1 152 


') 1 .745 


2.98 






- 




0.2730 


1.6582 


0.4072 


3.00 


0.10 


32. 


26. 




3.00 


19. 


0.058 


0.2559 


0.2885 


0.1 152 


')4.742 


2.98 


1.258 


10' 


2.620 


10" 


0.2767 


1.7177 


0.4283 


3.00 


0.30 


110. 


7fi. 






67. 


0.210 


0.2560 


0.3042 


0.1201 


75.96') 


2.97 


5.556 


10= 


2.764 


10 


0.2972 


2.6327 


0.5873 


3.00 


0.50 


185. 


135. 


'175 H87 


3.00 


126. 


0.394 


0.2571 


0.536') 


0.1 H80 


77.623 


2.97 


4.049 


10= 


2.891 


10" 


0.3054 


3.6'J71 


0.685 1 


3.00 


0.60 


206. 


166. 


339.749 


3.00 


160. 


0.497 


0.2585 


0.7fil3 




76.023 


2.97 


3.232 


10= 


2.627 






4.4546 


0.7298 


3.00 


0.70 


271. 


1 99. 


342.777 




198. 


0.607 


0.2605 


1 . 1 276 


0.3 1 53 




2.97 


3.342 


10= 


3.161 


10" 


0.3100 


5.3672 


0.7647 


3.00 


0.80 


295. 


236. 


345.579 


3.00 


247. 


0.751 


0.2618 


1.4397 


0.3635 


81.176 


2.97 


2.836 




3.(I2K 


10" 


0.3131 


6.6684 


0.8059 


3.00 
3.00 


0.90 
0.95 


389. 
436. 


281, 
307. 


347.185 
348.099 


3.00 
3.00 


278. 
272. 


0.830 
0,810 


0,2628 
0J2668 


1.6815 
2.4237 


0.3960 
0.4789 


77.906 
83,829 


2.97 
2,97 


2.043 
1.985 


'!)= 


2.637 
2.562 


10" 
10" 


0.3140 
0.3157 


8.4523 
9.5293 


0.8444 
0.8633 


250 


m 


3*4 


25 


440 910 


250 


19 


0059 


02559 


02885 


1152 


181 413 


248 


8 377 


10= 


1 .934 


10" 


02792 


1 6?m 


04025 


2.50 


0.30 


105. 


75. 


514.068 


2.50 


64. 


0.208 


0.2559 


0.2989 


0.1184 


159.866 


2.48 


3.598 


10= 


1.975 


10" 


0.3003 


2.5368 


0.5645 










slqfi^O 




Is^ 
















10= 


1.854 


10" 


ft 






9 


60 


198 


160 




9 SO 




408 


2582 


6761 


2205 




7 48 


1 739 


10= 


1.673 


10" 




4 1915 


6971 


250 
2.50 
2.50 


70 
0.80 
0.90 


236 
309. 

372. 


192 
228. 
271. 


555 058 
559.658 
563.501 


250 
2.50 
2.50 


191 
240. 
263. 


614 
0.764 
0.827 


2597 
0,2608 
0.2626 


9429 
1.1566 
1.4650 


2759 
0.3122 
0.3591 


155 846 
154.259 
154,899 


248 
2.48 
2.48 


1 744 
1.494 
1.299 


10= 
10= 
10= 


1.787 
1.787 
1.792 


10" 
10" 
10" 


3152 
0.3153 
0.3167 


5 1187 
6.0331 
7.3960 


7370 
0.7590 
0.7946 


2.50 


0.95 


418. 


298. 


565.360 


2.50 


265. 


0.831 


0.2663 


2.1052 


0.4356 


153.547 


2.48 


1,332 


10= 


1.914 


10* 


0.3190 


8.3689 


0.8175 


2.00 
2.00 


0.00 
0.10 


0. 

29. 


0. 
25. 


805.801 
839.531 


2.00 
2.00 


0. 
19. 


0.000 
0.061 


0.2559 
0.2559 


0.2885 
0.2885 


0.1152 
0.1152 




















2.00 


0.30 


88. 


73. 


947.948 


2.00 


63. 


0.214 


0.2560 


0.2951 


0.1172 




















2.00 


0.50 


150. 


127. 


1007.382 


2.00 


120. 


0.412 


0.2571 


0.4409 


0.1596 




















2.00 
2.00 


0.60 
0.70 


213. 
253. 


156. 
188. 


1022.614 
1034.130 


2.00 
2.00 


152. 
191. 


0.519 
0.647 


0.2584 
0.2600 


0.5977 
0.7980 


0.1985 
0.2414 




















2.00 


0.80 


309. 


224. 


1045.843 


2.00 


237. 


0.795 


0.2605 


0.9126 


0.2611 




















2.00 


0.90 


366. 


269. 


1(M8.349 


2.00 


243. 


0.805 


0.2645 


1 .4343 


0.3397 




















2.00 


0.95 


426. 


292. 


1054.389 


2.00 


246. 


0.817 


0.2671 


1.8486 


0.3842 




















1.70 
1.70 


0.00 
0.10 


0. 

27. 


0. 
25. 


1284.504 
1 344.255 


1.70 
1.70 


0. 
19. 


O.OOO 
0.063 


0.2559 
0.2559 


0.2885 
0.2885 


0.1152 
0.1152 




















1.70 


0.30 


83. 


71. 


1519.179 


1.70 


62. 


0.215 


0.2560 


0.2920 


0.1163 




















1.70 


0.50 


163. 


123. 


1616.978 


1.70 


117. 


0.413 


0.2573 


0.4062 


0.1494 




















1.70 
1.70 
1.70 
1.70 
1.70 


0.60 
0.70 
0.80 
0.90 
0.95 


199. 
239. 
287. 
342. 
371. 


152. 
183. 
219, 
262. 
284, 


1646.526 
1666,675 
1679,739 
1677,683 
1695,775 


1.70 
1.70 
1.70 
1.70 
1.70 


149. 
188. 
236. 
241. 
238. 


0.525 
0.655 
0.811 
0.817 
0.812 


0.2587 
0,2604 
0.2609 
0.2657 
0,2687 


0.5352 
0.6918 
0.7790 
1.2567 
1.5883 


0.1820 
0.2164 
0.2330 
0.3066 
0.3451 





















3.504 
3.122 

2.332 
2.475 
2.849 



13.90 
13.41 

13.50 
13.80 
13.76 



End of C- burning 











0.3214 


2,7453 


0.6739 


8.667 


10' 


7.305 


10' 


0.3318 


3.2222 


0.8165 


1.475 


10' 


2.689 


10^ 


0.3647 


5.3239 


1.0507 


6.487 


10" 


1.772 


10* 


0.3845 


8.8557 


1.2540 


5.534 




1.791 


10* 


0.3885 


1 1 .8407 


1.3410 


4.66H 




1.754 




0..'959 


16.')753 


1.5184 
1 X231 





10.8t4 










0.2977 


2.460t4 


7.630 


10.68 


8.494 


10^ 


1.264 


10^ 


0.3060 


2.8601 


6.173 


10.30 


2.397 


10^ 


6.628 


10* 


0.3293 


4.5785 


5.168 


11.34 


4.826 


10" 


1.910 


10* 


0.3499 


7.7164 


4.116 


11.39 


4.067 


10* 


1.837 


10' 


0,3579 


10.9273 


5.920 


11.24 


3.295 


10* 


1.931 


10' 


0,3569 


15.8098 


0.078 


11.25 


2.889 


10* 


2.812 


10' 


0.3592 


24.2797 


0.059 


10.89 


3.432 


JO" 


2.774 


10' 


0.4114 


56.7852 


0.040 


10.98 


3.258 


10" 


2.320 


10' 


0.4299 


77.6552 


0.001 


8.61 










0.2736 


1.8993 


0.046 


8.69 


1.589 




9.650 


10' 


0.2773 


2.3089 


0.053 


8.65 


5.740 




5.775 


10' 


0.3004 


3.7268 


0.047 


8.69 


3.279 




5.744 


10' 


0.3146 


6.6121 


0.008 


8.75 


2.615 




6.749 


10' 


0.3196 


9.5706 


0.060 


8.87 


1.91 1 




3.621 


10' 


0.3213 


13.8242 


0.043 


8.89 


1.509 




4.548 


10' 


0.3274 


22.0803 


0.044 


8.86 


1.367 


10" 


3.918 


10' 


0.3565 


46.2713 


0.057 


8.84 


1.436 


10" 


3.451 


10' 


0.3808 


67.9385 



2 0796 
1 .9853 
0.5905 
0.7254 
0.8818 
1.0572 
1.1534 
1,2025 
1.2414 
1.5159 
1.6540 



0.4749 
0.5711 
0.7339 
0.8869 
0.9622 
1.0176 
1.0901 
1.2773 
1.4335 



Notes. Same data as in Tab. 2 
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Georgy et al.: Populations of rotating stars 
Table 4. Main parameters of A-B stars at Z = 0.002 = Zsmc 



LnJ i)f H-biimiiig 



End of He-burning 



Mq 




"'km 


VMS 




M 

Mq 






I'^uii 


N/C 
mass fract. 


N/0 




M 
Mq 


P, 






I'surf 


N/C 
mass fract. 


N/0 


re 
kyr 


M 
Mo 


day 


a/a 




N/C 
mass fract. 


N/O 


15.00 
15.00 


0.00 
0.10 
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168. 


834.069 


2.00 


169. 


0.509 


0.2579 


1.2870 


0.3249 


198.499 


1.97 


1.964 


10* 


1.867 


10* 


0.3050 


6.0322 


0.7007 




















2.00 


0.70 


268. 


201. 


837.165 


2.00 


209. 


0.618 


0.2616 


1.9881 


0.4031 


0.000 


1.98 


7.316 


10* 


3.162 


10= 


0.3067 


8.0971 


0.7341 




















2.00 


0.80 


322. 


238. 


843.153 


2.00 


259. 


0.753 


0.2656 


2.9834 


0.4785 






































2.00 


0.90 


396. 


282. 


849.985 


2.00 


285. 


0.810 


0.2678 


3.7471 


(1.5247 






































2.00 


0.95 


409. 


308. 


871.563 


2.00 


298. 


0.842 


0.2753 


5.6638 


0.6399 






































1.70 
1.70 


0.00 
0.10 


0. 

29. 


0. 
27. 


1053.592 
1 101.952 


1.70 
1.70 


0. 
21. 


0.000 
0.064 


0.2509 
0.2509 


0.2885 
0.2885 


0.1152 
0.1152 






































1.70 


0.30 


104. 


77. 


1253.898 


1.70 


71. 


0.226 


0.2513 


0.3202 


0.1248 






































1.70 


0.50 


173. 


133. 


1313.265 


1.70 


134. 


0.421 


0.2552 


0-6981 


0.2181 






































1.70 
1.70 
1.70 
1.70 
1.70 


0.60 
0.70 
0.80 
0.90 
0.95 


213. 
251. 
304. 
360. 
420. 


164. 
196. 
232. 
267. 
298. 


1337.826 
1337.166 
1538.864 
1444.859 
1701.731 


1.70 
1.70 
1.70 
1.70 
1.70 


169. 
210. 
265. 
303. 
272. 


0.525 
0.636 
0.796 
0.854 
0.798 


0.2592 
0.2636 
0.2809 
0.2802 
0.3111 


1.1158 
1.6972 
4.0962 
5.2415 
12.7300 


0.2898 
0.3599 
0.5528 
0.5449 
0.9124 







































Notes. Same data as in Tab. 2 
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